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Abstract 

We consider a system-reservoir model where the reservoir is modulated by 
an external noise. Both the internal noise of the reservoir and the external 
noise are stationary, Gaussian and are characterized by arbitrary decaying 
correlation functions. Based on a relation between the dissipation of the 
system and the response function of the reservoir driven by external noise we 
derive the generalized Kramers' rate for this nonequilibrium open system. 
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I. INTRODUCTION 



More than half a century ago Kramers proposed a diffusion model for chemical reac- 
tions in terms of the theory of Brownian motion in phase space [Q]. Since then the model 
and several of its variants have been ubiquitious in many areas of physics, chemistry and 
biology for understanding the nature of activated processes in classical [0-0], quantum and 
semiclassical pUj-P^ systems, in general. These have become the subject of several reviews 
0,T^-rH| and monograph [r^| in the recent past. 

In the majority of these treatments one is essentially concerned with an equilibrium 
thermal bath at a finite temperature which stimulates the reaction co-ordinate to cross the 
activation energy barrier. The inherent noise of the medium is of internal origin. This im- 
plies that the dissipative force which the system experiences in course of its motion in the 
medium and the stochastic force acting on the system as a result of random impact from 
the constituents of the medium arise from a common mechanism. From a microscopic point 
of view the system-reservoir Hamiltonian description [p0|-p3| developed over the decades 
suggests that the coupling of the system and the reservoir co-ordinates determines both 
the noise and the dissipative terms in the Langevin equation describing the motion of the 
system. It is therefore not difficult to anticipate that these two entities get related through 
a fluctuation-dissipation relation ( these systems are sometimes classified as thermody- 
namically closed system in contrast to the systems driven by external noise in nonequilibrium 



statistical mechanics ||2^ ). However, when the reservoir is modulated by an external noise 
it is likely that it induces fiuctuations in the polarization of the reservoir. These fiuctuations 
in turn may drive the system in addition to usual internal noise of the reservoir. Since 
the polarization fiuctuations of the reservoir crucially depend on its response functions, one 
can further envisage a connection between the dissipation of the system and the response 
function of the reservoir due to the external noise from a microscopic point of view. 

In the present paper we explore this connection in the context of activated rate processes 
when the reservoir is modulated by an external noise. Specifically our object here is to 
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calculate the generalized Kramers' rate for the steady state of a nonequilibrium open system. 
Both the internal and the external noises are Gaussian, stationary and are characterized by 
arbitrary decaying correlation functions. While the internal noise of the reservoir is thermal, 
the external noise may be of thermal or non-thermal type. We consider the stochastic 
motion to be spatial diffusion limited and calculate the rate of escape in the intermediate 
to strong damping regime. We further mention that nonequilibrium, non-thermal systems 
have also been investigated phenomenologically by a number of workers in several other 
contexts, e.g., for examining the role of color noise in stationary probabilities [Q, properties 



of nonlinear systems nature of crossover p8[, rate of diffusion limited coagulation 



processes effect of monochromatic noise etc. While these treatments concern 
direct driving of the system by an external noise the present consideration is based on 
modulation of the bath. A number of different situations depicting the modulation of the 
bath by an external noise may be physically relevant. As, for example, we consider a simple 
unimolecular conversion (say, an isomerization reaction) from A — > B. The reaction is 
carried out in a photochemically active solvent under the influence of external fluctuating 
light intensity. Since the fluctuations in the light intensity result in the fluctuations in the 
polarization of the solvent molecules, the effective reaction field around the reactant system 
gets modified. Provided the required stationarity of this nonequilibrium open system is 
maintained (which is not difficult in view of the experiments performed in the studies of 
external noise- induced transitions in photochemical systems ||31[|) the dynamics of barrier 



crossing becomes amenable to the present theoretical analysis that follows. 

The remaining part of this paper is organized as follows: In Sec. II we discuss a system- 
reservoir model where the later is modulated by an external noise and establish an interesting 
connection between the dissipation of the system and the response function of the reservoir 
due to external noise. The stochastic motion in a linearized potential field is described in 
terms of a Fokker-Planck equation in Sec. III. Based on the traditional flux over population 
method we derive in Sec. IV the generalized expression for the Krmaers' rate of escape 



from a metastable well. The general theory is illustrated with a specific example in Sec. V. 



The paper is concluded in Sec. VI. 



II. THE SYSTEM-RESERVOIR MODEL : THE RESERVOIR MODULATED BY 

EXTERNAL NOISE 

We consider a classical particle of mass M is coupled to a heat bath of harmonic 
oscillators. The various modes of heat bath are perturbed by an external random force. The 
interaction between external force and bath co-ordinates is added as a new term {Hint) in 



the standard system-reservoir Hamiltonian of Zwanzig form ||2^. The total Hamiltonian is 
given by 

H=^ + V{x) + 1^1^ + m,u;Kq^ - a.{x)f^ + H^^t • (1) 

In Eq.(|T]), X and p are the co-ordinate and momentum of the system particle; {qi,Pi) are the 
variables associated with the i-th oscillator and uji and mi are the corresponding frequency 
and mass respectively. ai{x) measures the interaction between the particle and the bath. 
V{x) is the potential energy of the particle. Hint is assumed to be of the form 

1 ^ 

Hint = - ^Kiqie{t) . (2) 
^ i=i 

The coupling function Kj measures the strength of interaction and e(t) is the external noise 
which we assume to be stationary and Gaussian with zero mean, i.e., (e(t)) = 0. 

We now consider the interaction between the system and the heat bath to be linear, i.e., 

ai{x) = giX , (3) 



and then eliminate the bath degrees of freedom in the usual way p0| , p2| , |23| to obtain the 
following generalized Langevin equation 



X = V 

dV 

V 



- I dt' 7(t - t') v{t') + fit) + n{t) (4) 
while constructing Eq.(D we have set M and rrij equal to unity ] where 



N 

lit) =Y.9i COS Uit . (5) 

i=l 

fit) is the fluctuation generated through the couphng between the system and the heat bath 
and is given by 

N 

f{t)=^gi {[qi{0) - gix{0)] u'^ cosuJit + Vi{0) Ui sint^itj . (6) 

i=l 

f{t) thus depends on the couphng functions Qi and the normal mode initial conditions, 
which we assume to be canonically distributed. The statistical properties of f{t) can then 
be summarized in the following equations; 

(/(t)) = and (7a) 
{f{t)f{t'))=kBT^{t-t') . (7b) 

In Eq.(^), 7r(t) is a force term driven by the external noise e{t) and is given by 

n{t) = - l\{t-t') t{t') dt' , (8) 
Jo 

where 

N 

ip{t) =^giKi LOiSinUit . (9) 

i=l 

The statistical properties of 7r(t) are determined by the normal mode distribution of the 
bath, the coupling of the system with the bath, the coupling of the bath with the external 
noise and the external noise itself. It may also be noted that the external noise e{t) and the 
internal noise f{t) are independent of each other because they have been generated from 
different origin. Eq.(^ is reminiscent of the familiar linear relation between the polarization 
and the external field where vr and e play the role of the former and the later, respectively. 
ip{t) can then be interpreted as a response function of the reservoir due to external noise 
e{t). 

In the continuum limit 7(t) and (p{t) reduce to the following forms 

7(t) = j du V{uj) g^{uj) cos cut (10) 



and 

(p{t) = J did ^(iu) K(iu) uj g{uj) sin tut . (11) 
where T^iuj) is the density of modes of the heat bath. If the couphng functions 5'(to') and 



k{uj) are assumed to be of the following forms to obtain a finite result in the continuum 
limit, 

g{uj) = ^ and k{uj) = y/r^ u 

/TV U 



where 5^0 and kq are constants and Tc is the correlation time, i.e., ^ is the cutoff frequency 

of the harmonic oscillators then the expressions for 7(t) and ip{t) reduce to 

2 

= ^ f dujV{u) cosut (12) 

and 

(^{t) = go Kq J duj V{uj) UJ sin cut . (13) 
From the above two relations, we obtain 

? = -^W • (14) 

Eq.(|l4D is an important content of the present model. This expresses how the dissipative 
kernel 7(t) depends on the response function y9(t) of the medium due to external noise e(t) 
[ see Eq.(H) ]. Such a relation for the open system can be anticipated in view of the fact 
that both the dissipation and the response function crucially depend on the properties of the 
reservoir especially on its density of modes and its coupling to the system and the external 
noise source. In what follows we shall be concerned with the consequences of this relation 
in terms of the Langevin description in the next section ( Eq.(|l^) ). 

III. GENERALIZED FOKKER-PLANCK DESCRIPTION OF THE LINEARIZED 

MOTION 



We now consider the system to be a harmonically bound particle of unit mass and of 
frequency ujq. Then because of Eq.(|) the Langevin equation (§) becomes 
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V = —UnX 



dt' -f{t - t') v{t') + f{t) - / dt' ip{t - t') e(t') 



(15) 



The Laplace transform of Eq. (|T3|) allows us to write a formal solution for the displacement 
of the form 



x{t) = {x{t)) + / dt' hit - t') fit') - —T,tol I dt' hit - t') eif) 



90 



— T, I dt' /l2(t - t') tit') 

9o Jo 



(16) 



where we have made use of the relation ([T^) explicitly. 
Here 

{xit)) = Xxit)xiO) + hit)viO) 



(17) 



with a;(0) and f (0) being the initial position and initial velocity of the oscillator, respectively, 
which are nonrandom and 



l-ool hir) dr 



(18) 



The kernel hit) is the Laplace inversion of 



his) 



s'^ + 7(,s) s + uJq 

where, 7(5) = /q°° e~** 7(t) dt, is the Laplace transform of the friction kernel 7(t), and 

d^hit) 



(19) 



h2it) 



dt^ 



(20) 



The time derivative of Eq. (|T^) yields 



vit) = Mt)) + /* dt' hit - t') fit') - —t,ujI f dt' hit - t') tit') 
Jo go Jo 

rt 



^r, / dt' hit~t') tit') 

go Jo 



(21) 



where 



(vit)) = -ujlhit)+viO)hiit) 



(22) 



and 

Next we calculate the variances. Prom the formal solution of x{t) and v{t), the explicit 
expressions for the variances are obtained which are given below; 

<^iit) = (Kt) - {x{t))r) 

= 2 tdh h{h) f' dt2 h(t2) {f(h)f(t2)) 
Jo Jo 

+2 l—T,iun f dh h{h) r dt2 h{t2) (e(ti)e(i2)> 

\go J Jo Jo 

+2 I'^tA [ dh h2{ti) [ ' dt2 h2{t2) {e{ti)e{t2)) 
\ Qo I Jo Jo 

f dh h{h) r dt2 h2{t2) {e{h)e{t2)) , (25) 
JO Jo 



and 



,^0 



"lit) = {[<•(') - {vmi 

ft rti 



2 / dh hih) I ' dt2 hit2) {fitl)f{t2)) 

Jo Jo 

+2 (—T^COn dh hih) dt2 h{t2) {e{ti)€{t2)) 



\9o 

+2 I—tA t dti h^{ti) r dt2 h^{t2) {e{h)e{t2)) 
\go J Jo Jo 
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+2(— rj [ dh hiih) [ ' dt2 h3{t2) {e{h)e{t2)) , (26) 
\go I Jo Jo 



'^Ut) ^ Mt) - {x{tmv{t) - {vit))]) 

- l^lit) (27) 

where we have assumed that the noises f{t) and e{t) are symmetric with respect to the time 
argument and have made use the fact that f{t) and e{t) are uncorrelated. 



8 



Due to the Gaussian property of the noises f{t) and e(t) and the hnearity of the Langevin 
equation (p!5D, we see that the joint probabihty density p{x,v,t) of the oscillator must be 
Gaussian. The joint characteristic function associated with the density is 



pifi, p, t) = exp \ i{x{t))fi + i{v{t))p - i + 2al^{t)pfi + cTl^it)p^ 



(28) 



Using the method of characteristic function p^|35| and the above expression (^) we find the 
general Fokker-Planck equation associated with the probability density function p{x, v, t) for 



the process (p!5|) 



where 



(29) 



7(t) 
Tit) 



dt 



InT(t) 



-h{t) hi{t) + hl{t) 



T{t) 



and 



hit) 



cUq / dr /i(r) 

JO 



+ hit) . 



The functions 0(t) and tpit) are defined by 



(pit) = uj^it)(Tl^ + -fal^ + , 



vv 



(30) 



where the covariances are to be calculated for a particular given noise process. 

For the internal noise processes Carmeli and Nitzan and Hanggi [^] have shown 
that for several models the various time dependent parameters ouoit), 7(t), etc. do exist 
asymptotically as t —>■ oo. The above consideration shows that hit), hiit), etc. do not 
depend on the nature of the noise but depend only on the relaxation 7(t). 

We now discuss the asymptotic properties of 0(t) and ipit), which in turn are dependent 
on the variances cr^^(t) and o"^^(t), as t oo since they play a significant role in our further 
analysis that follows. 
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From Eqs. (^) and (pG]), we may write 



2(e)/ 



and 



where 'i' denotes the part corresponding to internal noise f{t) and 'e' corresponds to the 
external noise e{t). Since the average velocity of the oscillator is zero as t ^ oo we see from 
Eq.(p2D that h{t) and hi{t) must be zero as t ^ oo. Also from Eq.(p!7D we observe that 
the function Xx{i) must decay to zero for long times. Hence, from Eq.(p!8D we see that the 
stationary value of the integral of h{t) is I/uJq, i.e., 



POO \ 

/ hit) dt = ^ . 

Jo LUq 

Now, crl^^^t) and cr'^lj\t) of Eqs.p5|) and (pBf ) can be written in the form 

Jo Jo 



(31) 



2 dr h{T) - h\t) - iul dr /i(r)' 



(32) 



and 



ali'\t) = 2 1 dh I dt2 hiih) (/(ti)/(t2)) 



(33) 



From the above two expressions [ Eqs.(p2D and ( p3| ) ] we see that 



and aJ!^Hoo) = keT . 



tug 



(34) 



It is important to note that these stationary values are not related to the intensity and 
correlation time of the internal noise. 

We next consider the parts, crli^\t) and cr^^'^''(t), due to the presence of the external noise. 
The Laplace transform of Eq. ([iBD yields the expression 
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x{s) - (x(s)) = h{s)f{s) -TcUjlh{s)e{s) -Tcs'^h{s)e{s) (35) 

where 

= |1 _ ^2^1 3;(o) + /i(s)t;(0) . (36) 

From the above equation (^) we can calculate the variance o"^^ in the Laplace-transformed 
space which can be identified as the Laplace transform of Eg. (P3|) . Thus, for the part cr^^^^(t) 
we observe that, o'^'^xK^) contains terms like (^^TcUJQh{s)^ (e^(s)). Since, we have assumed 
the stationarity of the noise e(t), we conclude that if C(0) exists [ where C{t — t') = (e(t)e(t')) 
], then the stationary value of cr^(f)(t) exists and becomes a constant that depends on the 
correlation time and the strength of the noise. Similar argument is also valid for cr'^l^^t). 
Summarizing the above discussions we note that, 

(i) the internal noise-driven parts of cr^.^.(t) and (T^„(t), i.e., a^^?^ and a^^*^ approach the fixed 
values which are independent of the noise correlation and the intensity as t cxd, 

(ii) the external noise driven parts of variances also approach the constant values at the 
stationary {t oo) limit which are dependent on the strength and the correlation time of 
the noise. 

Hence we conclude, following the Refs. [36,37] and our preceding discussions that even 
in presence of an external noise the coefficients of the Fokker-Planck equation (ppj ) do ex- 
ist asymptotically and we write its steady state version for the asymptotic values of the 
parameters as, 

dp n dp d d'^p d'^p . . 

- + + ^9^^"^^ + "^^"^^d^ + ^^^^ 9^ = ° ' (^^^ 

where, uJq, 7, (p{oo), ip^oo), etc. are to be calculated from the general definition (|30| ) for the 

steady state. 

The general steady state solution of the above equation (^) is 



Pstix,v) = -^exp 



2 — 2 

V UJqX 



2Do 2(Do + ^(oo)) 
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(38) 



where 



fl„ = ^ (39) 

7 



and Z is the normahzation constant. The solution ( p8| ) can be verified by direct substitution. 
The distribution (|38| ) is not an equihbrium distribution. This stationary distribution for the 
nonequihbrium open system plays the role of an equilibrium distribution of the closed system 
which may, however, be recovered in the absence of external noise term. 

Some further pertinent points regarding rate theory for nonequihbrium systems may be 
in order. It is well-known that the equilibrium state of a closed thermodynamic system with 
homogenous boundary conditions is time- independent. The open, e.g., the driven system 
on the other hand, may show up the complicated spatio-temporal structures or may settle 



down to multiple steady states . The external noise may then induce transitions between 
them. It is, however, important to realize that these features originate only when one takes 
into account of the nonlinearity of the system in full and the external noise drives the system 
directly. Secondly, in most open nonequihbrium systems the lack of detail balance symmetry 
gives rise to severe problem in the determination of the stationary probability distribution 
for multidimensional problems [^. At this juncture three points are to be noted. First, for 
the present problem we have made use of the linearization of the potential at the bottom 
and at top of the barrier (as is done in Kramers' and in most of the post-Kramers' 
development |]7|-p|, p!6| , ^ ) which precludes the existence of a multiple steady states. Second, 



the external noise considered here drives the bath rather than the system directly. Third, 
the problem is one-dimensional. Thus, an unique stationary probability density which is an 
essential requirement for the mean first passage time or fiux over population method (as in 
the present case) for the calculation of rate and which is readily obtainable in the case of 
closed equilibrium system, can also be obtained for this open nonequihbrium system. 
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IV. KRAMERS' ESCAPE RATE 



We now turn to the problem of decay of a metastable state. To this end we consider as 
usual a 'Brownian particle' moving in a one- dimensional double well potential V{x) ||T],|TB[. 
In Kramers approach the particle coordinate x corresponds to the reaction coordinate 
and its values at the minima of the potential V{x) denotes the reactant and product states. 
The maxima of V{x) at Xf, separating these states corresponds to the activated complex. All 
the remaining degrees of freedom of both reactant and solvent molecules constitute a heat 
bath at temperature T. Our object is to calculate the essential modification of Kramers' 
rate when the bath modes are perturbed by an external random force when the system has 
attained a steady state. 

Linearizing the motion around barrier top at the Langevin equation (^) can be 

written down as 



y = v 



v = ujly- dt' 7(t - t') v{t') + fit) + 7c{t) , (40) 





where, y = x — Xb and the barrier frequency uj^ is defined by 

V{y) = H - \u;h' ; ^6 > . (41) 

Correspondingly the motion of the particle is governed by the Fokker-Planck equation ( p9D 

where, the suffix '6' indicates that all the coefficients are to be calculated using the general 
definition (^) for the barrier top region. 

It is apparent from Eqs. (|37D and ( ^21) that since the dynamics is non-Markovian and 
the system is thermodynamically open one has to deal with the renormalized frequencies 
u!q and uJb near the bottom or top of the well, respectively. Following Kramers |1| we make 
the ansatz that the nonequilibrium, steady state probability pb, generating a nonvanishing 
diffusion current j, across the barrier is given by 
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Pb[X,V) 



exp 



' V{x) 



(43) 



where 



lb 



(44) 



V{x) is the renormahzed hnear potential as 

1 

2 
1 



V{x) = V{xo) + ^ijl{x-Xof 



V{x) = V{xb) - -uJbix - XbY 



near the bottom 
near the top 



(45) 



with uJq, ojI > 0. The unknown function ^{x,v) obeys the natural boundary condition that 
for X ^ oo, C,{x,v) vanishes. 

The ansatz of the form ( ^31) denoting the steady state distribution is motivated by the 
local analysis near the bottom and the top of the barrier in the Kramers' sense For 
a stationary nonequilibrium system, on the other hand, the relative population of the two 
regions, in general, depends on the global properties of the potential leading to an additional 
factor in the rate expression. Although because of the Kramers' type ansatz which is valid 
for the local analysis, such a consideration is outside the scope of the present treatment, we 
point out a distinctive feature in the ansatz (|43|) compared to Kramers' ansatz. While in 
the latter case one considers a complete factorization of the equilibrium part (Boltzmann) 
and the dynamical part, the ansatz ( ^31) incorporates the additional dynamical contribution 
through dissipation and strength of the noise into the exponential part. This modification of 
Kramers' ansatz (by dynamics) is due to nonequilibrium nature of the system. Thus unlike 
Kramers', the exponential factors in ( ^3]) and in the stationary distribution (|38| ) which serves 
as a boundary condition are markedly different. While a global analysis may even modify 
the standard Kramers' result our aim here is to understand the modification of the rate due 
to modulation of the bath driven by an external noise, within the perview of Kramers' type 
ansatz. The internal consistency of the treatment, however, can be checked by recovering 
the Kramers' result when the external noise is switched off. 
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From equation (^21), using ( ^3]) we obtain the equation for $,{y,v) in the steady state in 
the neighborhood of Xb, the equation 



Moo) ] 
Db y dy 



Db + Moo 



^ + 0,(oo)p + ^,(oo)£|- = o . (46) 



dv 



We then make use of the following transformation 



u = V + a y , y = X — Xb 



dvdy 



where a is a constant to be determined. We obtain from Eq.(^) 



{0fe(oo) + aipbioo) 



Db 



Db + Moo, 



-ujIx 



MpoY 

Db 



du 



(47) 



Putting 



Db 



-2 , - , h , Moo) 



-Xu 



(48) 



Db + Moo)"" ' ' V ' Db _ 
(where A being another constant to be determined) we obtain the ordinary differential equa- 
tion for ^(m) 



where 



A 



d'e. di 

-— + Au— = 

au^ du 



A 



0fe(oo) + aMoo) 

and the two constants A and a must satisfy the simultaneous relations 

Db 



— Xa 



Db + ^^6(00 



-A = 76 + a 1 + 



MooY 



Db 



This implies that, the constant a must satisfy the quadratic equation 

Db 



Db + Moo) 2 , - 
a + 7b a 



Db + Moo] 



(49) 



(50) 



(51a) 
(51b) 



(52) 



which allows 
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a± 



The general solution of Eq. (^9|) is 



{-7fc±\/¥+4^^} • 



(53) 



^{u) = F2 I exp 
/o 



dz + Fi 



(54) 



where Fi and F2 are constants of integration. We look for a solution which vanishes for 
large x. For this to happen the integral in should remain finite for |m| +00. This 
implies that A > so that only a_ becomes relevant. Then the requirement pb{x,v) —>■ 
for X —> +00 yields 

^ (55) 



F^ — F2 \ / — r 

^' 2A 



Thus we have 



7r 



2A+yo 



and correspondingly 

Pb{x,v) = F2 




2A + io 



A^^ 



exp 



+ 



V{x) 



2Db Db + Moo) 



(56) 



The current across the barrier associated with this distribution is given by 



+00 



V pb{x = Xb, v) dv 



which may be evaluated using (K3) and the linearized version of V{x)^ namely, V{x) 



V{xb) - hujiix - XbY, as 



j = F2 



2ti 



A + D 



\ 1/2 
— Db exp 



V{xb) 



Db + ^b(oo) 



(57) 



To determine the remaining constant F2 we proceed as follows. We first note that as x —00 
the pre-exponential factor in Eq. (|56D reduces to the following form 



FJ.... 



A 



(58) 



We then obtain the reduced distribution function in x as 
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Pb{x -oo) = 2nF2 ( — 1 exp 



V{x) 



where we have used the definition for the reduced distribution as 



(59) 



p{x, v) dv 



Similarly we derive the reduced distribution in the left well around x ~ using Eq.(p8| 
where the linearized potential is V{x) = V{xo) + jUq^x — xqY, 



Pst[x) = — \/2nDo e oo+^i°°) exp 



with the normalization constant 1/Z 



UJ'^{X - XqY 



g-Do+i/'(oo) 



Z 27r^Z}o(/^o + ^(oo)) 
The comparison of the distributions ( p9D and (|60|) near x = xq, i.e., 



Pst[Xo) = Pb[Xo) 



(60) 



(61) 



gives 



DbJ 271 J271 {Do + (oo)) 



exp 



V{xo) 



Hence from 



Db + ipb{oo)_ 

, the normalized current or the barrier crossing rate k is given by 



(62) 



Db 



A 



1/2 



exp 



Db + ipb[oo) 



(63) 



27T {Do + ^{oo)y/^ \l + ADb. 

where Eq is the activation energy, Eq = V{xb) — V{xQ). Since the temperature due to internal 
thermal noise, the strength of the external noise and the damping constant are buried in the 
parameters Dq, Db, ipo, ipo and A the generalized expression look somewhat cumbersome. 
We point out that the subscripts '0' and '6' in D and ip refer to the well or barrier top region, 
respectively. Eq. (|63D is the central result of this chapter. The dependence of the rate on the 
parameters can be exposed explicitly once we consider the limiting cases. 
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V. EXAMPLE : (5-CORRELATED EXTERNAL NOISE AND 
ORNSTEIN-UHLENBECK INTERNAL NOISE 



We consider a particular case where the external noise e{t) is (5-correlated and the internal 
noise is an Ornstein-Uhlenbeck (0-U) process, i.e., 

(e(t)e(t')) = 2D6{t-t') 

and 

{fit) fit')) = glksT 



2, 



One can trace the origin of the above correlation function for internal noise by considering 
a Lorentzian type frequency distribution of the normal mode variables and Eq.(^. Conse- 
quently, from the fluctuation-dissipation relation we derive the dissipative kernel as 

-|t-t'|/Tc 

lit-^)=gl- . (64) 

It should be noted that for Tc — 0, the above noise processes become 5-correlated. 

The correlation functions as given above and the dissipative kernel provide the required 
quantities for the calculation of Pst(x, v) \ Eq.(|38|) ] and the generalized rate expression (|63|). 
Thus we have [see the Appendix] 

(/,(oo) = 7(A;bT + Dk^) and ?/'(oo) = . (65) 

which gives 

Z^o = ^ = W + DkI . (66) 

7 

Hence from equation ( pSf ) we see that the steady state distribution is given by 



Pstix,v) = ^exp 



(67) 



2ikBT + DK'o) 

since for a Markovian process lDq = uJq. From Eq.(p7D we see that, the steady state prob- 
ability density in our example does not depend on the correlation time of noise but does 



depend on the strength and couphng of the external noise. The result is in agreement with 
the one obtained earlier by Bravo et. al. 

We now return to our generalized rate expression Eg. (|63D . For the present case we have 
computed [see the Appendix] 







Da 



Ijj 



b ' 



giikBT + D^l) 



! 



and 



(68a) 
(68b) 
(68c) 
(68d) 

(68e) 



Using all these values, we obtain from Eq.(|63|) 



k 



27lUJh 



1/2 



9l 



En 



exp 



ksT + DkI 



(69) 



If we put the external noise intensity D equal to zero, i.e., when external noise is absent, 
the above expression reduces to usual Kramers rate expression with = 7. We note here 
that Dn^/kB defines a new effective temperature characteristic of the steady state of the 
nonequilibrium open system. As expected this temperature is function of the strength of 
external noise intensity D and the coupling of the external noise to the bath modes. 



VI. CONCLUSIONS 

Based on a system-reservoir microscopic model where the reservoir is modulated by an 
external, stationary and Gaussian noise with arbitrary decaying correlation function, we 
have generalized the Kramers' theory to calculate the steady state rate of escape from a 
metastable well. The main conclusions of this study are as follows; 

(i) We have shown that since the reservoir is driven by the external noise and the dissipative 
properties of the system depend on the reservoir, a simple connection between the dissipation 
and the response function of the medium due to the external noise can be established. 
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(ii) This connection is important for realising the stationary state of the thermodynamically 
open system characterized by an effective temperature of the reservoir, which depends on 
the strength of the external noise. 

(iii) Provided the long time limit of the moments for the stochastic processes pertaining to 
the external and internal noises characterized by arbitrary decaying correlation functions 
exist, the expression for generalized Kramers' rate of barrier crossing for the open system 
we derive here is fairly general. The expression assumes simple forms in the specific limiting 
cases. 

The creation of a typical nonequilibrium open situation by modulating a bath with the 
help of an external noise is not an uncommon phenomenon in applications and industrial 
processing. The external agency generating noise does work on the bath by stirring, pump- 
ing, agitating, etc., to which the system dissipates internally. In the present treatment we 
are concerned with a nonequilibrium steady state characterized by an effective temperature 
which signifies a constant throughput of energy in contrast to thermal equilibrium defined 
by an constant temperature. We believe that these considerations are likely to be important 
in other related issues in nonequilibrium open systems. 

ACKNOWLEDGMENTS 

SKB in indebted to Council of Scientific and Industrial Research (CSIR), Government 
of India for financial support. 



APPENDIX A: CALCULATION OF VARIANCES 



The Laplace transform of 7 in Eq.(|64D is given by 




and subsequently, we have 
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rcS^ + s^ + {LjiTc + gi)s + iui 

For Tc 7^ 0, the above expression can be simplified to 

~ s + a 

^ ^ s^ + as^ + bs + Co 

where 

a = — , = uJq -\ and cq — — . 

''"c ''"c ''"c 

We define the following characteristic quantity 

Q = ^ \ + 

108 27 27 6 4 
and distinguish the cases : Q > 0, Q = and Q < 0, for which we have different forms of 
h{t). For the case Q > [ we do not give here the explicit expressions for the cases Q = 
and Q < ] we find that the inverse Laplace transform of h{s) reads 

h{t) = cie-^i* + C2e-^2* sm{pt + a) (Al) 

where the coefficients Ci, C2, Ai, A2, (3 and a are given by 

A,^-A-B + ^ , (A2a) 
A, = ^{A + B) + ^ , (A2b) 

f3^^{A-B) , (A2c) 

= 2A2-A,-c^ ' ^^2d) 

-^-"^'^^"g^-A,^^"''^ , (A2e) 

a6 Co ^ 



C2 = -"j-Kd - A2)2 + fiY'^ and (A2h) 
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Here we note that for a physically allowed solution Ai, A2 must be positive. Since by Eq.(|T9D 
h{t) depends on the memory kernel 7(t) which is of decaying type and all the moments, in 
general, reach asymptotic constancy as shown in Sec. Ill, these quantities are positive (which 
depends on the correlation time Tc, the strength of the noise and other potential parameters) 
which may be checked (after some algebra) by considering the limiting cases such as Tc ^ 
and Tc — large. 

Substituting Eq. (|55|) into the expressions for variances [external noise is 5-correlated] , 
namely into (p5|) and (|2^) we have after some lengthy algebra 

where 



a: 



2(0 



it) = kBT{c2R+^^ 



2c2e-^2t 



2 - 2uj^qC2R 

1 - UJ'^C2R + 



2 D 2a;^ci 



+ e 



Lolci 



2 



Ai 



Ai V ^1 , 

[A2 sin(/5t + a) + l3cos{(3t + a)] 



-2ciC2e-('^^+^^)* sin(/?t + a) 
A2/?a;2cie-2^^* 



2, ,2^2^--2A2t 



(A2 + /52)2 
0.2(2/3^ 



+ 



■sin2(/3t + a) - 



(Ai + /32)2 



Ai) 



(Ai + /32)2 



- 1 



sin^(/5t + a) 



(A3) 



with 



Ai + /?2 



(A2 sin a + /? cos a) 



(A4) 



and 



2D (^^r,j [ciiuj', + At + 2ujlA\)lA{t) 

+cl{ujt, + (A^ - (3^ - A^^^l + 2ujI{AI - (3^)}lB{t) 

+2c,C2{ujt + A?(A^ - + ^2(^2 ^ ^2 _ (32^^^I^it) 

-2cl[3A2{Al -f3^ + ujl)lD{t) + AclP^AllE{t) - 4ciC2/?A2(A? + ujD^m 



(A5) 
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Here the /'s are defined by 



J^(t) = / dt , (A6a) 

Jo 

lB(t)= I e-^'^^'sm^((3t + a) dt , (A6b) 
Jo 

Ic(t)= f e-(^^+'^^)* sin(/?t + a) rft , (A6c) 
Jo 

lD(t)= I e-^^'^ sin 2(pt + a) dt , (A6d) 
Jo 

Ij^{t) = [ dt and (A6e) 

Jo 

lF{t)= r e-(^^+^^)* cos(/5t + a) tit . (A6f) 



Similarly 



<^lit) = cr^it) + cr^it) 



where 



-/3A2C^e-2^2* sin 2(/3t + a) + (A^ + 00^- (3^)cle-^'^^' sin^if^t + a) 
+e-(^^+^^)*{2ciC2(u;o + ^1^2) sm{pt + a) - 2Ai/5ciC2 cos(/?t + a)}] (AT) 



and 



\t)=2D [c\Al{ul + AlfUit) 

+cl{{ul + A^ - W'f^l - {^l + 3A^ - 
+2C1C2A1A2K' + A?)(^o - 3/3' + A2)/^(t) 

-C2/?A2(3A2 - + ^2)(^2 _ 3^2 ^ ^2)j^(^) 

+c^/32K'-/32 + 3A^)J^(t) 

-2ciC2/?Ai(A? + uI){ujI + 3A2 - /32)J^(t)] (A8) 



where, J's are defined in Eg. ( |A6a[ - [A61| ). The explicit expression for o"^^(t) can be derived 



from Eg. (P?]) . From Egs.(0) and (pDj), we calculate the stationary values of the variances 
as follows, 
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where 



and 



XX 



aZX' (oo) = ksT C2R + -r— 



Cl 

Ai 



2-a;2(c2i?+^ 



+ ^27^(r^ + 4/3^^2) 



(A9) 



4(Ai + /32, 

A? + ^2 



{/3(r2 - 4/32 - 4Alf3) sm2a - A2{r^ + 4/3T - AAjp^) cos2q;} 



+2CiC2 



(Ai + A2)2+/52 



X ({r(Ai + A2) + 2/?2A2} sin a + /3{r - 2A2(Ai + A2)} cos a)] 



(AlO) 



where 



(All) 



Again substituting the values of ci, C2, Ai, A2, a, /3, F and we obtain after a lengthy 
calculation the following result 



(A12) 



Similarly 



+c^^{/3^(^o' + 3A^ - 13^ + A^((^o^ + A^ - 



4(Ai + /32) 

X ({/?A2(u;o^ + A2 - 3/32) (a;2 + 5^2 + ^2^ ^ ^3(^2 ^ 3^2 _ ^2^2^ 2a 

+ {A2/52(o^2 ^ 3^2 _ ^2)^^2 _ ^2 _ 5^2) ^ ^3^^2 ^ ^2 _ 3^2)2| ^^^^a) 

Ai(A^+a;g) 

+^'^'^Ai + A2)2 + /?2 

X ({A2(Ai + A2)(o;2 + A2 - 3f3^) - ^^u'^ + SAj - P^)}sma 



+{PA2{u;l + Al- 3/?2) - p{Ai + A2){u;^o + 3A^ - /?')} cos a)] . (A13) 
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Again putting the values of all the parameters we have, 

al{^)^kBT + D4 . (A14) 

Clearly, 

alioo) = . (A15) 

The variances al^{oo), (j^^(oo) and (T^^(oo) yield 0(oo) and ip{oo) and other relevant quan- 
tities. 
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